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Abstract 

^ . 

• , We study the optimal control of a steady-state dead oil isotherm prob- 

nI . lem. The problem is described by a system of nonlinear partial differen- 

tial equations resulting from the traditional modelling of oil engineering 
within the framework of mechanics of a continuous medium. Existence 
and regularity results of the optimal control are proved, as well as neces- 
sary optimality conditions. 
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1 Introduction 

We are interested in the optimal control of the steady-state dead oil isotherm 
problem: 

-/S.tf{u) = div {g{u)Vp) in fl, 

- div id{u)Vp) = / inn, 



\on = 0- 



(1) 

.Plao no- 
where n is an open bounded domain in R^ with a sufficiently smooth boundary. 
Equations ([Ij serve as a model for an incompressible biphasic flow in a porous 
medium, with applications to the industry of exploitation of hydrocarbons. The 
reduced saturation of oil is denoted by u, and p is the global pressure. To 
understand the optimal control problem that we consider here, some words 
about the recovery of hydrocarbons are in order. For a more detailed discussion 
about the physical justification of equations ((T|) the reader is referred to [H [TBI 
[It] and references therein. At the time of the first run of a layer, the flow of 
the crude oil towards the surface is due to the energy stored in the gases under 
pressure in the natural hydraulic system. To mitigate the consecutive decline of 
production and the decomposition of the site, water injections are carried out, 
well before the normal exhaustion of the layer. The water is injected through 
wells with high pressure, by pumps specially drilled to this end. The pumps 
allow the displacement of the crude oil towards the wells of production. The 
wells must be judiciously distributed, which gives rise to a difficult problem of 
optimal control: how to choose the best installation sites of the production wells? 
This is precisely the question we deal within this work. These requirements lead 
us to the following objective functional: 

Jiu,pJ)^^\\u-U\\l + ^\\p-P\\l + ^\\f\\ll, (2) 

where 2 > go > 1 and /3i > is a coefficient of penalization. The first two terms 
in ^ make possible to minimize the difference between the reduced saturation of 
oil u, the global pressure p and the given data U and P, respectively. Our main 
goal is to present a method to carry out the optimal control of ([!]) with respect 
to all the important parameters arising in the process. More precisely, we seek 
necessary conditions for the admissible parameters u, p and / to minimize the 
functional J. 

Theoretical analysis of the time-dependent dead oil problem with different 
types of boundary and initial conditions has received a significant amount of 
attention. See [3] for existence of weak solutions to systems related to ([T]), 
uniqueness and related regularity results in different settings with various as- 
sumptions on the data. So far, optimal control of a parabolic-elliptic dead oil 
system is studied in [IBj . Optimal control of a discrete dead oil model is consid- 
ered in |18| . Here we are interested to obtain necessary optimality conditions 
for the steady-state case. This is, to the best of our knowledge, an important 
open question. 



Several techniques for deriving optimality conditions are available in the 
literature of optimal control systems governed by partial differential equations 
[51 [SI [ini [m [131 [m • in this work we apply the Lagrangian approach used with 
success by Bodart, Boureau and Touzani for an optimal control problem of the 
induction heating [2], and by Lee and Shilkin for the thermistor problem [7]. 

The motivation for our work is threefold. Firstly, the vast majority of the 
existing literature on dead oil systems deal with the parabolic-elliptic system. 
Considering that the relaxation time for the saturation of oil u is very small, 
the time derivative with respect to the saturation is dropped. Hence we get 
the system ([T|). Such a steady-state dead oil model represents a reasonably 
realistic situation where we neglect the time derivative. Secondly, some technical 
difficulties when dealing with system ([T|) arise and rely on the fact that there is 
no information on the time derivative of the reduced saturation of oil as well as 
on the pressure. As a result, one cannot use directly the standard compactness 
results to obtain strong convergence of sequences of solutions in appropriate 
spaces. This is in contrast with J15j . where a fully parabolic system is considered. 
Thirdly, the choice of the cost function ^ for this time dependent problem 
seems to be quite appropriate from the point of view of practical applications. 

The paper is organized as follows. In Section [2] we set up notation and 
hypotheses. Additionally, we recall two lemmas needed in the sequel. Our 
main results are stated and proved in the next two sections. Under adequate 
assumptions {HI) and {H2) on the data of the problem, existence and regularity 
of the optimal control are proved in Section [3l In Section [H making use of 
the Lagrangian approach and assuming further the hypothesis (ff3), we derive 
necessary optimality conditions for a triple (u, p, /) to minimize ([2|) among all 
functions {u,p,f) verifying ^. We end with Section [5] of conclusion. 

2 Preliminaries 

The following assumptions are needed throughout the paper. Let g and d be 
real valued C^-functions and (p he a. C^ function. It is required that 

(HI) < ci < d{r), g{r), ip{r) < ca; C3 < d'(r), ip'{r), ip"{r) < a for all r G M, 
where Ci, i — 1, . . . , 4, are positive constants. 

(H2) U, P e L^{fl), where U, P -.n^R. 

(H3) \ip"'{r)\ <cfor allreR. 

Henceforth we use the standard notation for Sobolev spaces: we denote 

II • Up = II • IIlp(o) for each p e [1, 00] and 

W^ = W^{n) := {u e LP{^), Vu e LP{n)} , 
endowed with the norm ||u||vi/i(n) = ll"llp + ll^^llpJ 

w^ = w^{n) := {u e w^m, W'u G LP{n)} , 



with the norm Ilu||jy2/Qi = ll'"llvyi(o) + ll'^^^ll • Moreover, we set 

V :- wl{n)- 

T := {! e L^'^in)] , 

ol 

where Wp (^) is the interior of Wp{V,). 

In the sequel we use the following two lemmas in order to get regularity of 
solutions. 

Lemma 2.1 f |12[|20j ). Let fl C i?" be a bounded domain with a smooth bound- 
ary. Assume that g € (L^(r2))" and a G C(Jl) with miuj^ a > 0. Let u be the 
weak solution to the following problem: 

—V • (aVw) = V • 5 in £7 
u = on dfl. 

Then, for each p > 2, there exists a positive constant c* , depending only on n, 
n, a and p, such that if g G {L'^{fl))", then 

\\^u\\p<c*i\\g\\p+\\yu\U). 

Lemma 2.2 ([6,). For any function u e C"(fi)n W2 {^) n W2 (^) ^'*e'~e exist 
numbers Nq and go such that for any g < go there is a finite covering of O 
by sets of the type flg{xi), Xi G 0., such that the total number of intersections 
of different 0,2g{xi) — VLC\ B2g{xi) does not increase Nq- Hence, we have the 
estimate 

4 / „IL,I|2 2a /^||V72„,||2 , 1 ||v7„,l|2 



||VW||: < C h||^„(f,) g'" [\\y'u\\^ + ^ IIV7.II2 

3 Existence and Regularity of Optimal Solutions 

In this section we prove existence and regularity of the optimal control under 
assumptions (HI) and {H2) on the data of the problem. 

3.1 Existence of Optimal Solution 

The following existence theorem is proved using Young's inequality together 
with the theorem of Lebesgue and some compactness arguments of Lions [8]. 
The existence follows from the fact that J is lower semicontinuous with respect 
to the weak convergence. Recall that along the text constants c are generic, and 
may change at each occurrence. 



Theorem 3.1. Under the hypotheses (HI) and (H2) there exists a q > 1, 
depending on the data of the problem, such that the problem of minimizing ([2]) 
subject to (d} has an optimal solution [u,p,f) satisfying 

Proof Let (m™,p'", /"*) G W^2^(^) x ^ x L'^'i°{Q) be a sequence minimizing 
J{u,p, /). Then we have that (Z™) is bounded in L'^''°{i}). By the second equa- 
tion of ([T|) governed by the global pressure and a general result of elliptic PDEs 
[1], under our hypotheses we have that Vp™ is bounded in L^'?(J7). Writing now 
the first equation of ([!]) as 

-div(^'(u™)V3(w'")) = div((7(M'")Vp™) 

and using Lemma [^TTl we obtain Vm™ G L'^''{n). Hypotheses allow us to express 
again the first equation of ([T]) as 

-^'(w'")Au™ - ¥>"(«") I Vw^p = div(g(w'")Vp™). 

Hence, 

h"'lk2(n) <c, 

where all the constants c are independent of m. Using the Lebesgue theorem and 
compactness arguments of Lions [S], we can extract subsequences, still denoted 
by (p"), (m") and (/"), such that 

u'"" ^u weakly in Wg{n), 

p" ^ p weakly in W^qi^), 
/™^7weaklyinL2'Jo(^)_ 

Then, by Rellich's theorem, we have 

p™ -> p strongly in L^{il). 

Therefore, by using these facts and passing to the limit in problem ([T]), it follows 
from the weak lower semicontinuity of J with respect to the weak convergence, 
that the infimum is achieved at [u, p, f). □ 

3.2 Regularity of Solutions 

Regularity of solutions given by Theorem 13.21 is obtained using Young's and 
Holder's inequalities, the Gronwall lemma, the De Giorgi-Nash-Ladyzhenskaya- 
Uraltseva theorem, an estimate from [S], and some technical lemmas that can 
be found in [B]. 



Theorem 3.2. Let (u,p, /) be an optimal solution to the problem of minimizing 
([2]) subject to (H)). Suppose that (HI) and (H2) are satisfied. Then, there exist 
a > such that the following regularity conditions hold: 

u,peC^{n), (3) 

u,peWl{n), (4) 

u,pew^{n), (5) 

ueCim, (6) 

u&wl^{n), pewig^in), (7) 

where qo appears in the cost function ([2]) • 

Proof. Firstly, ([3]) is an immediate application of the general results of [SI [51 HH] ■ 
To continue the proof of Theorem 13. 2i we need to estimate ||Vu||4 in function 
of ||Vp||4. Taking into account the first equation of ([1]), it is well known that 
u G Wlifl) (see [5]) and 

llVulU < c||Vp||4. (8) 



Using Lemma |2.2[ we have, for any g < go, that 



||Vp||t<c|b||^„(^)g2"|||Vp||^ + l||Vp| 



Therefore, we get (|3]) for an eligible choice of g. Using ([5]), we obtain that 
u S Wl{i^). On the other hand, by the first equation of ([Ij and the regularity 
^, we have that u e W|(ri). Moreover, it follows, by the fact that u € W|(0), 
that p e Wi{n). Using again (g]) and the fact that W2(f7) -^■ Ci(n.), the 
regularity estimate (|B]) follows. Finally, the right-hand side of the first equation 
of (HD belongs to L-* (17) ^ L^'>° (n) as 2qo < 4. Thus, by © weget m e VF|^^(rj). 
Since / G i^'^°(J7), the same estimate follows from the second equation of the 
system ([T]) for p. D 

4 Necessary Optimality Conditions 

We define the following nonlinear operator corresponding to ([T]): 

F -.W xW xT — > H 
{u,p,f)^F(u,p,f) = Q, 

where 

p(„ „ f\- f -^v(w) - div(.g(u)Vp) 
r[u,p,j)-^ - div (d(u)Vp) - / 

Due to the estimate 

lkllH.i^(0)<c||«IU2^(n), ^veW^^in), l<q<2 

2-q 



(see [S]), hypothesis (HI) and regularity results (Theorem I3.2p . we have 

(y9'(u)Au, ip"{u) I Vw|^ , g'{u)\7uVp, d{u)VuVp G L^ (f7) C L^'^in). 
Thus, it follows that F is well defined. 

4.1 Gateaux Differentiability 

Theorem 4.1. Lei assumptions (HI) through (H3) hold. Then, the operator F 
is Gateaux differ entiahle and its derivative is given by 

SF{u,p, /)(e, w, h) = —F {u + se,p + sw, f + sh) 1^^^ = {SFuSF2) 

- div {ip'{u)'Ve) — div {(p"{u)e'Vu) — div (g(u)Vw) — div {g'(u)eVp) 
- div (rf(u) Vw) - div {d'{u)eVp) - /i 

/or aZ/ {e,w,h) ^W xW y. T. Furthermore, for any optimal solution (u,p,jj 
of the problem of minimizing ^ among all the functions {u,p,f) satisfying ([I}, 
the image of 5F (u,p,f) is equal to H. 



To prove Theorem 14.11 we make use of the following lemma. 

Lemma 4.2. The operator 6F{u,p,f) : W x W x T — > H is linear and 
bounded. 

Proof. We have for all {e,w,h) ^W xW xT that 

6pF2{u,p, f){e,w, h) = — div {d{u)\7w) — div ((i'(u)eVp) — h 
= -d{u)Aw - d'{u)Vu ■ Vw - d'{u)eAp - d'{u)Ve ■ Vu - d'{u)eVu -Vp-h, 

where 5pF is the Gateaux derivative of F with respect to p. Then, using hy- 
pothesis (HI), we obtain that 

\\SpF2{u,p, /)(e, w, h)\\2g < ||Vw||2, + c\\Aw\\2g 
+ c\\Vu ■ Vw\\2g + c\\eAp\\2q + c||Ve • Vu||2, + c\\eVu ■ Vphg + \\h\\2q. (9) 

We proceed to estimate the term ||eVw- Vp||2ij. Similar arguments can be applied 
to the remaining terms of ([9|). We have 

\\e\/u-\/p\\2q < ||e||oo||VM-Vp||2g 

< l|e||oc||VM||^||Vp||4 

2-q 

< c||M||vK|b||w||e||n/. 
Then, 

\\SpF2{u,pJ)ie,w,h)\\2g <ci\\u\\w, Iblk, II/IIt) i\\e\\w + Mw + ||/^||t) • 

(10) 



On the other hand, 

6uFi{u,p,f){e,w,h) 

= — div {(p'{u)Ve) — div ((^"(u)eVu) — div {g{u)\'w) — div {g'{u)eVp) 
= ~ip'{u)Ae - Lp"{u)Vu ■ Ve - ip"{u)eAu - ip"{u)\/e ■ Vu - ip"'{u)e\\/u\^ 
— g{u)Aw — g'{u)\7u ■ Vw — g'{u)eAp — g'{u)\/e ■ \7p — g" {u)e\Iu ■ Vp, 

where 5uF is the Gateaux derivative of F with respect to u. The same arguments 
as above give that 

\\5uF^{u,pJ){e,w,h)h, < cdl^llH^, Ibllv^, II/IIt) (||e|k + \\M\w + WHt) ■ 

(11) 
Hence, by ([TO]) and (ITT]) . 

\\5F{u,pJ){e,W,h)\\Hy,Hy.T <c(\\u\\w. \\p\\w . ||/||t) (||e|k + \\w\\w + ll/illx) . 

Consequently the operator SuFi{u^p, f) is hnear and bounded. D 

Proof of Theorem \4- 1\ In order to show that the image of 6F{u,p, f) is equal 
to H, we need to prove that there exists a. {e,w,h) G W x W x T such that 

— div ((^'(u)Ve) — div (ip" {u)eVu) — div {g(ll)Vw) — div {g' {u)eyp) — a, 
- div (d(u)VM;) - div (d'(M)eVp) - /i = ^, 

for any (a, /3) G H. Writing the system (fT2|) for h = as 

—Lp'{u)Ae — 2(p" (u)Vu ■ Ve — cp" {u)eAu — (p'" {u)e\\7u\ 
—g{u)Aw — g'{u)Vu ■ Vw — g'{u)eAp — g'{u)'S/p ■ Ve — g"{u)e\/u ■ Vp — a, 
-d{u)Aw - d'(u)Vu ■ Vw - d'{u)eAp - d'{u)Vu ■ Ve - d'{u)eVu -Vp^ P, 

elan = 0> 
= 0, 

(13) 



w|an = 



it follows from the regularity of the optimal solution (Theorem 13. 2p that 

(/?"(u)Au, (^"'(u)|Vup, g'{u)Ap, g"{u)Vu ■ Vp, d'{u)Ap, d'{u)Vu ■ Vp e L^'">{n), 
(/3"(m)Vu, 5'(u)Vu, g'{u)Vp, d'{u)Vu e L*®(f}). 

By general results of elliptic PDEs [6l [H [19] , there exists a unique solution of 
system ^^ and hence there exists a (e,w,0) verifying p^ . We conclude that 
the image of (5i^ is equal to H. D 



4.2 Optimality Condition 

We consider the cost functional J: WxWxT^M.^ and the Lagrangian C 
defined by 

'C(u,p, /,pi,ei) = J{u,pJ) + (f{u,pJ), f ^^ 

where the bracket (•, •) denote the duahty between H and H' . 

Theorem 4.3. Under hypotheses (H1)~(H3), if (u,p, f) is an optimal solu- 
tion to the problem of minimizing ([2]) subject to ^, then there exist functions 
(eijPi) G Vl^2 (f2) X 14^1 (rj) satisfying the following conditions: 

div ((/j'(u)Vei) - d'{u)Vp ■ VpT - ^" (u) Vu • VeT - g'{u)Vp ■ VeT = m - C/, 

div (d(u)Vpr) + div {g{u)Ve{) =p- P, 

Pi\dn = 0, 

2go/3i|7|'*-7-pr. (14) 

Proof. Let (u, p, / j be an optimal solution to the problem of minimizing ([2]) 
subject to ([TJ. It is well known (cf., e.g., [5]) that there exist Lagrange multipli- 
ers {pi,ei) e H' verifying 6(^u.pj)^ (u,p, /,pr, el) {e,w, h) — for all {e,w, h) G 
W xW xT, with 5(u,pj)t^ the Gateaux derivative of £ with respect to (u,p, /). 
We then obtain 

((u - U)e + {p- P)w + 2go/3i |7|'*"'7/i) dx 
n 



+ / (— div ((^'(u)Ve) — div ((^"(u)eVu) — div {g{u)\Iw) — div {g'(u)eVp)) ei da; 

+ / (- div (d{u)Ww) - div {d'{u)eVp) - h)pldx = 
Jo 

for all {e,w,h) G W x W x T. This last system is equivalent to 



{{u — U)e — div {d'{u)e'Vp)pi — div ((y9'(u)Ve) ei 

— div ((/3"(u)eVu) el — div {g' {u)eVp) el) dx 



+ / ((p — -P)w — div ((i(u)Vw)pi — div ((7(u)Vu)) ei) dx 

+ / (2(?o/3i|7P'°-'7/i -pT/i) d2; = 



for all {e,w,h) € W x W x T. In other words, 

[ {{u^U) + d'{u)\7p-VpT 
Jn 

- div ((/3'(u)VeT) + ip" (u) Vu ■ VeT + g'{u)\/p ■ Ve^)edx 

+ / {{p- P) -div (d(u)VpT)~ div {g{u)WeT))wdx 
Jn 

+ f {2qol3i\f\^''°-^Jh-pTh) dx = (15) 
Jn. 

for all {e,w,h) ^ W x W x T. Consider now the system 

div ((^'(u)Vei) - d'{u)Vp ■ Vpi - ip"{u)\7u ■ Vei - g {u)S/p ■ Vei ^u-U, 
div (d(u)Vpi) + div {g{u)Vei) ^p- P, 

eilan = Pilan = 0- 

(16) 
It follows again, by [6l[8l[T9], that ((T6)) has a unique solution (ei,pi) G ^^2^(51) x 
W|(i7). Since the problem of finding {e,w) & W x W satisfying 

— div (<y9'(u)Ve) — div {ip" {u)eVu) — div {g{u)\'w) — div {g'{u)e'Vp) 

— sign(ei — eT) — div {d{u)\/w) — div {d'{u)eVp) (17) 

= sign {pi - pT) 

is uniquely solvable on W2q x W2q, choosing h = in ([T5|) . multiplying ([TBI) by 
(e,w), integrating by parts, and making the difference with (jlSp . we obtain 

(- div {ip'{u)Ve) - div {ip"{u)eVu) - div (.g(u)Vw) - div (.g'(u)eVp)) 
^^ (18) 



X (ei — ei) dx + (- div ((i(u)Vw) - div (d'(u)eVp)) (pi — pi) da; = 
Jn 

for all (e, w) S M^ x W^. Choosing (e, w) in (jTS)) as the solution of system (jT?)) . 
we have 

/ sign(ei - eT)(ei - eT)(ix+ / sign{pi -pT){pi -pi)dx ^ 0. 
Jn Jn 

It follows that 61=6]" and pi = pT. On the other hand, choosing (e, w) = (0, 0) 
in ((T5|) , it follows (ITil) , which concludes the proof of Theorem 14.31 D 

5 Conclusion 

In this paper, we considered the optimal control of a steady-state dead oil 
isotherm problem with Dirichlet boundary conditions, which is obtained from 
the standard parabolic-elliptic system, where the relaxation time for the reduced 
saturation of oil is very small. The main purpose was to prove existence and 
regularity of the optimal control and then necessary optimality conditions. The 
proposed method is based on the Lagrangian approach. 
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